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ABSTRACT 

An adaptive control scheme for aircraft was studied to find the 
cause of limit cycles which occurred under large input commands and to 
find a method for eliminating the oscillations. The complexity of the 
system equations dictated that all analytical studies be performed on 
simplified versions of the adaptive control scheme. After exhaustive 
analysis of the cause of limit cycles in the simplified systems and an 
investigation of possible fixes for this problem, it was decided that 
the oscillations were an inherent feature of the control scheme 
resulting from the servo and actuator lags. 

A second adaptive control scheme was thus developed which 
obtained more feedback information downstream of the aircraft actuator 
and series servo. A comparison of the two adaptive systems was made 
under similar flight conditions for the F-4 aircraft and the modified 
adaptive controller was shown to be practical and free of limit cycle 
oscillations by analog simulation. 
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I. INTRODUCTION 


During the past few years, a new control concept has developed in the • 
field of automatic control. The concept of self-adaptive control made 
possible the capability of inherently satisfying desired system perform¬ 
ance criteria in the presence of a changing environment. An adaptive 
control system can change its parameters through an internal process of 
measurement and evaluation of dynamic performance. 

Although self-adaptive control systems were found to be successful 
in most respects, they do exhibit undesirable characteristics. Con¬ 
trollers which detect limit cycles, either intentionally, as in the 
Honeywell control system for the X-15, or unintentionally, as in the 
General Electric control system for the F-lll, can be upset by air 
turbulence unless special precautions are taken which reduce the 
effectiveness and increase the complexity of the system. One such con¬ 
trol system, a nonvarying-C* adaptive control scheme for aircraft was 
proposed by Rang, Reference 1. The desired system performance criterion, 
the Boeing C*-Criterion, specifies that the time response of the quantity 
C* fall within a prescribed envelope, for all speeds and altitudes, thus 
the name nonvarying-C* controller. This control system, as other similar 
adaptive systems, exhibits undesirable limit cycles under high input com¬ 
mands. An analog study was conducted in an effort to find the cause of 
these limit cycles and provide a method or methods of eliminating the 
oscillations. The study was accomplished at the Naval Postgraduate 
School Hybrid Computer Facility, Monterey, California, during the period 
November, 1968, through June, 1969. 
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The hybrid computer facility consists of a high speed digital com¬ 
puter model 9300 manufactured by Scientific Data Systems of Santa Monica, 
California, and an electronic analog computer model Cl 5000 manufactured 
by COMCOR, Ind., of Anaheim, California, together with extensive soft¬ 
ware. Various displays are available including a X-Y recorder, a six 
channel recorder, and two general purpose high performance input/output 
consoles which have cathode ray tubes and light pens as well as keyboards. 

The 9300-Computer is a high speed, general purpose digital computer. 
The various input/output devices include: teletype, paper tape, magnetic 
tape, light pen, and analog interface. The system is operated directly 
by the user and utilizes FORTRAN IV language. 

The Cl 5000 is a large scale, general purpose electronic analog/ 
hybrid computer. It is equipped with summers, integrators, inverters, 
diode function generators, electronic resolvers, quarter square multi¬ 
pliers, D/A converters, A/D trunks, comparaters, relays, manual and servo 
potentiometers. 

The adaptive control system simulation used only the analog computer; 
however, for future study on such problems the hybrid capability of the 
computer could be utilized by programming the digital machine to auto¬ 
matically set the potentiometers for each new flight condition or for 
different aircraft under study. 
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II. THE NONVARYING-C* ADAPTIVE CONTROLLER 


The original nonvarying-C* adaptive controller was developed beginning 

with the short period aircraft perturbation equations of motion. 

e=M0 + Mct + M-a + M.6 
q a a o e 

n = U (6 - a) = - Z a - Z.6 
o a 6 e 

where the aircraft coefficients M . M , M-, NL, Z , and 1 £ are constants 

q a a o or o 

representative of a specific flight condition. It is desirable to pro¬ 
vide uniform response in a changing environment without using measure¬ 
ments of altitude, dynamic pressure, or angle of attack. Further, angle 
of attack is difficult to measure and contributions to angle of attack 
provided by gusts are difficult to distinguish from the motion of the 
aircraft; therefore, angle of attack is eliminated from the equations. 


e = M e + M 
q a 


n 


+ M- 

0 - fl- 

Z 

Z e 

a 

U 

a 

a 

o_ 


+ M.6 

o e 


n = 


a 


9_U o 


- Z .6 . 

o e 


If the following substitutions are made. 


Z*M 

M - — 
6 Z 


a 


B = M + M- 

q q a 

M M- 

e n ‘ Z " U 
a o 


the fundamental equation in the study is derived as 

« ' V + + 

The desired system performance criterion, the Boeing C*- criterion, is 


the sum of the normal force at the pilot's seat, n + £p0, and a constant 


15 















multiple of pitch rate, V c 9, thus 

C* = n + JL0 + Vj. 

P c 

The normal acceleration term dominates at high speeds and the pitch rate 
term dominates at low speeds. If C* is required to be a constant mul¬ 
tiple of the input command, it may be seen that the C*~criterion is in 
exactly the same form as the fundamental equation, 

v c • i k 

6 = - -jp- 6 - 4— n - -jr— c 

*p x p p 

0 = i q e + $ n n + e fi 5 e 

where coefficients B q > B p , and 6^ are negative for all flight conditions. 
V_ 


Now set 


V' * p 




k_ 

i 


and the C*-criterion becomes 

0 = B„6 + B„n + $»C. 
q n 6 

If a control function could now be found whose gains vary so that for each 

different flight condition the control system behaved as if the funda- 

• • ** . “ — 

mental aircraft equation 0 - B q 0 - $ n n - B fi 6 e = 0 were replaced by the 
system performance criterion 0 - B q 0 - B n n - B^C = 0, then the system 
would be adaptive. If we neglect the actuator, servo, and actuator model 
in the analysis for the present, this control function may be taken as 


and 


f = r 0 + r n + r.C 
q n o 

5 = C + f. 

e 


The configuration is shown in Figure 1. 
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An error signal detects the amount the system criterion deviates 


from zero. 


e = 6 - 6 q e - 6 n n - 3 g C. 


Substituting the fundamental aircraft equation and the control function 
yields the following error relationship upon which the derivation of the 
variable gains is based; 



Thus the following ideal values of the variable gains will yield zero 
error: _ 



It is necessary to derive equations for these gains which will drive 
the system error to zero. Shipley, Reference 2, suggests a qradient 
technique which has become standard in such derivations. The calcula- 


tions are: 2v - B q + 6 { r q - 8 q 2 + jr 8„ + 6 { r n - B n 


2 



If the variable gains are chosen to make the qradient dv/dt < 0, then as 
time proceeds, v and hence e approach zero. 
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dr 

Choose: I, - = K 0G 

dt q 

^=KnG 
dt n 

dT 

~5T = K S CG 


Thus 


3T ' V G 


Observing that B fi is always negative and setting G equal to e or sgne, 
the gradient dv/dt is consistently negative. Thus at each flight condi¬ 
tion, the error will approach zero and the variable gains will approach 
their ideal values. 

Due to the complex nature of the equations of this adaptive control 
scheme, simplified versions of the scheme were developed which could be 
exactly and completely studied analytically. As a first cut, the 
simplest system was considered. The fundamental equation was reduced to 
the following: 

e' = e s s,. 


The entire system should now behave like the model or criterion, 0 = 0^C. 
The error equation reduces to e = 0 - B fi C and the control function 
reduces to f = r^C. Since 8^ = C + f (neglecting servo), the relation¬ 
ship for the ideal value of the variable gain, r^, can be calculated, 


e = 0 - BgC 

S ' Ve " S 6 C 
c - 6 S (1 + r 6 )C - 6 { C 




for e 5 0 


as was the case for the complete system. 
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Therefore, the controller is stable if 2?w > - (B g < 0). 

Therefore, for a given servo, c and u are fixed and the system may be 
driven unstable by sufficiently increasing K g or C Q . Since C Q is a 
squared term, the system stability is largely determined by this quantity. 
It is now obvious why the limit cycles are most predominant under large 
input commands. 
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III. ANALOG SIMULATION FOR THE C*-CONTROL SYSTEM 


The Cl-5000 analog computer was used to simulate the simplified con¬ 
troller. The second order servo was modeled by the circuit shown in 
Figure 3, where C-j and C2 were decade capacitor boxes. The values of 
z, and to for the servo could be provided by appropriately setting C-j and 
C 2 from the log-log graph in Figure 3. The simplified controller was 
first operated with G=e, which provided a linear system. The system was 
run without the servo, C-j and C 2 equal zero,, and the convergence of the 
gain to its ideal value was verified, Figure 4. The decade capacitor 
boxes were then set to give neutral stability, or a simulated limit 
cycle. Figure 5. The controller was altered by setting G = sgne and 
replacing a multiplier by a relay on the analog computer, see Figure 6. 
With the same settings on the decade capacitor boxes or any more unstable 
settings, this controller would exhibit limit cycles as shown in Figures 
7 and 8, instead of a divergent instability. 

A second simplified version of the nonvarying-C* controller was then 
studied to note the increase in complexity toward the full adaptive 
system and to verify the existence of limit cycles in an unstable non¬ 
linear system. The equations which describe this system now include the 
0 terms, but assume B g = 3^. The equations are more complex and include 
more of the terms of the complete controller, therefore yielding results 
which better approximate the complete control system. The set of equa¬ 
tions are: 

5 - V + V. 
e ■ « - V ‘ V 

« = c + f = c + r e 


23 







































































































































































































































fa 

fa 










f 



























e 

Rir/ 1 n na 




-- 























— 

— 







_ 


















_ 

















j 

-j 

















— 

J 




















_ 

i 









_ 





——j 




















£ 









— 

ft 

m 
















~ 

— 

Si 








B 

id 




















r* 

u 









& 

V 























1 A1 









±Ui a 

1 


-9* 




0 

2v/ line 


106 

2v/£i 


ne 


-it- 


e 

2 v/line 


FIGURE k. ANALOG TRACE FOR r g C- SYSTEM, WITHOUT SERVO C^O C 2 =0 
































































































































































































































































FIGURE 5. ANALOG TRACE FOR r\.C-SYSTEM, NEUTRAL STABILITY OSCILLATIONS 

d C 1 =.k7 C 2 =.0039 





























































































































































































































































































































FIGURE 6. ANALOG DIAGRAM FOR r C-SYSTEM 


a 













































































































































































































































































































































ro 

CO 


1 1 

_L___. 

\ 

i 


l i i 






r n! 

I 




! 










j 



















i 


































































c 

f, 









is 

v 

































































1 ,6 

2v/line 













1 

i 

i 

—i 

r 










FIGURE 8. ANALOG TRACE FOR r\C-SYSTEM, LIMIT CYCLE OSCILLATIONS (RELAY) 

C l=.47 C 2 =.0039 














































































































































































dF n 

dt 3 " = K q® e (G=E ^ 

5 + 2cu<$ + w 2 6 = uj 2 6 . 

e e e e 

The block diagram is shown in Figure 9. Combining the equations gives: 

dr . i- _ . i 

3r = K q 6 L (e q ' e 6 )e + (6 e ' C)B 6j 

6 e + 2?w6 e + w 2 6 e = w 2 (C + r q e) 

5 = i q 0 + 6 s 6 e , 

and applying the small perturbation theory 


r_ = r + r 

^o 


q '9n 'q 


5=6 +5 

e e 0 

e = e 0 + e 

yields the characteristic equation of the simplified r q 9 system: 
X 4 - (3 q - 2cto)X 3 - (2&o3 q - w 2 )X 2 

- (e q u 2 + e 6 w 2 r q ^ + 3 6 K q e o V)x + K q a) 2 e 0 2 e 6 e q = o , 

where r = A ~ —3. an d q = ■ 0 5 . 

3 0 B 6 0 B q 

Therefore the system is stable if 

- (6 y * - ^ C ° - ’ B i-)(2;i^ - «*)(£ - 2 W ) > 

q 3 q q 


K u) 2 C 2 3 p 3 - K w 2 C 2 3 3 

(V' * - 3 - ' ^ ~ ') + (B q - **“>*< q B V > 

and limit cycles will not occur. It is obvious from the preceding 
stability relationship for the r q 8-controller that the equations for 
the complete control system are too complex to analyze. The r q e-system 
is, however, sufficiently similar to the complete controller that most 
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FIGURE 9. BLOCK DIAGRAM FOR SIMPLIFIED r q 9- c C> WTR OL SYSTEM 




































of the analysis is applicable to both systems, as may be seen by 
analog simulation. 

The r^e-system was first operated without the servo and as the error 
approached zero, the variable gain approached its ideal value. Figure 10. 
The second order servo was again simulated by the circuit shown in Figure 
3 and with G = e, the decade capacitor boxes were set to give the system 
neutral stability. The resulting motion is shown in Figure 11. With 
G = sgne, a relay was used and the analog diagram for the system is 
shown in Figure 12. With sufficiently large inputs the system becomes 
unstable. This results in the limit-cycle oscillations shown in Figure 
13. 

The complete nonvarying-C* adaptive control system was then patched 
on the analog computer according to the diagram shown in Figure 14. In 
Reference 1, Rang shows that the equations for the controller are some¬ 
what more complicated than those derived in Chapter II resulting from 
the introduction of the actuator and actuator model. The complications 
resulting from the actuator and the second order servo will be shown in 
the following analog graphs. Figure 15 shows the free aircraft, that is 
the aircraft without the control function feedback loop, with no actuator. 
Figure 16 shows the free aircraft with the actuator. The effect of the 
actuator, considered as a first order lag 1/.1S+1, may be seen by com¬ 
paring 6 and 6 on the two preceding figures. Figures 17 and 18 show 
6 

aircraft parameters of the complete adaptive control system without the 
actuator or actuator model and Figures 19 and 20 show the system para¬ 
meters with the actuator and actuator model. As can be seen from these 
figures the variable gains do not approach the same value for the con¬ 
troller with and without the actuator. Thus the introduction of the 
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FIGURE 12. ANALOG DIAGRAM FOR f 0-SYSTEM 
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FIGURE 13. ANALOG TRACE FOR T 9-SYSTEM, LIMIT CYCLE OSCILLATIONS (RELAY) 

q C =1.28 C 2 =.005 
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FIGURE 14. ANALOG DIAGRAM FOR NONVARYING^C CONTROL SYSTEM 
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FIGURE 17. ANALOG TRACE FOR C -CONTROLLER, WITHOUT ACTUATOR 0^=0 C 2 =0 
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FIGURE 19. ANALOG TRACE FOR C -CONTROLLER, WITH ACTUATOR C^O C 2 =0 
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actuator has so complicated the control system that it no longer behaves 
as if its governing equation were exactly that of the criterion as the 
system error approaches zero. As soon as the second order servo is 
introduced into the system, the combined complications produce limit 
cycle oscillations as shown in Figure 21. It can thus be seen that the 
servo and actuator lags are complicating the system and causing the 
limit cycles. 

The results obtained for the complete nonvarying-C* adaptive con¬ 
troller without the actuator or servo may be verified by calculating 
values of the variable gains based on the aircraft parameters for the 
F-4 aircraft. Figure 18 shows the complete adaptive controller without 
the complications introduced by the actuator or the servo at a flight 
condition of Mach 2.0 at 50,000 feet. Flight Condition 3 on Table I 
contains aircraft parameters for the F-4 aircraft at this same flight 
condition. Therefore, using the values of the ideal variable gains 
calculated from these aircraft parameters and listed on Table II, a 
comparison may be made to the values obtained by analog simulation. 
Figure 18. The values correspond closely which indicates that at this 
flight condition the variable gains are approaching their ideal values 
and the error is approaching zero. Similar results were obtained for 
all flight conditions which indicates that the basic controller is 
operating properly and limit cycles will not appear. Only with the 
actuator and series servo in the system, which is the case in the actual 
aircraft control system, will limit cycles appear. 
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TABLE I 


AIRCRAFT AND CONTROLLER PARAMETERS, F-4 AIRCRAFT 


(units in radians, feet, seconds) 


FI ight 

Condition 

l 

2 

3 

Altitude 

Sea Level 

Sea Level 

50,000 

Mach No. 

0.200 

0.900 

2.000 

100/U 

0 

0.451 

0.100 

0.052 

-M /TOO 
a 

0.013 

0.380 

0.424 

-M-/10 

a 

0.026 

0.139 

0.004 

-M q /10 

0.043 

0.261 

0.048 

-M 6 /T 00 

0.028 

0.584 

0.149 

-Z /T 000 
a 

0.083 

2.420 

0.817 

-z 6 /iooo 

0.0131 

0.314 

0.102 

-e 6 /ioo 

0.200 

0.200 

0.200 

-3 q /l° 

0.300 

0.300 

0.300 

-103 n 

0.200 

0.200 

0.200 

K 6 /100 

0.999 

0.999 

0.999 

104K n 

0.500 

0.500 

0.500 

K q /10 

0.999 

0.999 

0.999 

A 

0.0146 

0.0143 

0.0519 

A 

0.690 

4.000 

0.520 

a 

2.580 

53.500 

9.610 

£ 

0.020 

0.020 

0.020 

0 
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TABLE II 

IDEAL VALUE OF VARIABLE GAINS 


Flight 

Condition 

1 

2 

3 

A1titude 

Sea Level 

Sea Level 

50,000 

Mach No. 

0.200 

0.900 

2.000 


Original Nonvarying 

-C* Controller 


V 

6.750 

-0.630 

1.080 

V 

0.900 

-0.019 

0.260 

* 

r 

n 

0.0021 

0.0001 

-0.0033 


Modified Nonvarying' 

-C* Controller 


V 

7.750 

0.370 

2.180 

V 

0.1155 

-0.050 

0.124 

4* 

r 

n 

0.00027 

0.00029 

-0.00159 
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IV. STUDY OF LIMIT CYCLES AND POSSIBLE FIXES 


In order to find fixes for the problem of limit cycle oscillations 
in a control system, the limit cycles must first be studied to find 
their amplitude and frequency. The combined equations for the control¬ 
ler using G = sgne are 
dr. 

cRT = K 5 Csgn t B ($ 6 e " 3 <S C ^ 

6 + 2ccoc$ + oj 2 6 = oj 2 (C + T.C) 

e b e o 

for the simplified r^C-system. Since the amplitude and frequency of the 
parasitic oscillations on the steady state are desired, apply the small 
disturbance theory by setting 



and 6 = <5 + 6 

e e e 

o 

Setting C = C Q and substituting in these perturbation equations gives 
the system's full variational equations 


dr. 
_ 0_ 

dt 

d 2 6 

dt 2 " 


K.C 
o 0 


+ 


sgn [6 { « a ] 
d6 

e , >£ 


aj 2 r .c „ 
0 0 


which may be combined to yield 

d 3 6 d 2 6 d6 _ A 

—^ + 2?w —^ + u 2 / - w 2 C 2 K.sgn [B.6 ] = 0. 

dt 3 dt 2 dt o 6 3 6 e 

The method of Harmonic Balance, Reference 3, allows higher order 
harmonics to be disregarded if it is assumed that <$ e = A sinftt + higher 
order harmonics. Making the proper substitutions for derivatives of 
<$ e in the preceding equation yields 
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AQ(Q 2 -w 2 )cosftt + 2cwAft 2 sinQt 
+ u 2 C 0 2 K g sgn[B 5 Asinnt] = 0. 

If the Describing Function for a relay with a deadband of 2 e q could 
be calculated, then it would be possible to eliminate the nonlinear term 
sgn[Asinftt]. According to Gibson, Reference 4, the Describing Function 
for such a relay is ksin (fit + 4>), where K = cos (sin -1 ) and 


A$, 


HA 

$ = 0 for this system. Making this substitution gives 

Aft(ft 2 - w 2 )cosftt + 

(2?toAft 2 - aj 2 C 0 2 K g ^ cos (sin -1 jg~)) sinftt = 0. 

Setting the coefficients of cosnt and sinfit equal to zero yields equa¬ 
tions for frequency and amplitude of the limit cycle oscillations> fi=aj 


A = 


|-[I+ (K 2 - 4e Q 2 /3/)*] 




where & = 2C 0 2 K s /ca>n. Figure 8 shows limit cycle oscillations for the 
r^C-system with C-j = .47 and C£ = .0039 which corresponds to w = 24 


rad/sec. and c = .14. With C Q = 1, 


.7, e = .2, and 3* = -10 


values of ft and A may be calculated from the preceding equations to be 
3.82 cycles/second and .1311 volts respectively. Figure 8 confirms 
these values on the trace for 106^. Since the frequency of the limit 
cycle oscillations is now known, it is known what frequencies may be 
filtered or where to apply a lead or lag compensation to possibly 
eliminate the limit cycles. 

An investigation of why the oscillations are present in the first 
place, results in a better understanding of the problem at hand. From 
the set of equations for the r^C- system in Chapter II, it was observed 


that the ideal values for 0 and r x were 0 = 3-.C and IV 

6 6 6 


= - i. if 

*5 


a oerturbed system is considered where 0 and r x are at some small 

o 
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disturbance quantity from their ideal values, the equations are 


c = c 


l 6 

*6 


r 5 = ^ - l + y 


£ = 9 - Vo 

and the small disturbance terms e and y are driving the system until 


e = 0, y = 0; therefore 0 = B X C and r r = r 

0 0 A 


Hence, from the system 


equations, 


8 . 


8gC Q + e = 8* [C 0 + (■=— 1 + y)C n ] and thus e = 8*C n y. The 


(To 1 


0 3 6 °- 

block diagram for the perturbed r^C-system is contained on Figure 22. 

This figure shows that there will be a maximum phase shift of -270 
degrees. It can also be shown that limit cycle oscillations occur at 
-180 degrees phase shift by the following simple diagram. 

e(t) 




Q 


nonlinear 

element 


linear 

element 


C(t) 


Assume, 


e(t) = Asintot 
C(t) = Aksin(wt + <J>) 

where k and $ are the amplitude and phase angle of the combined non¬ 
linear and linear elements. In order for a periodic oscillation to 
occur, it can be seen that e(t) must equal -C(t). Therefore, -Asin wt 
= Aksin(wt + $) and thus k = 1 and $ = -180 degrees, or the amplitude 
ratio must equal one and the phase shift must equal -180 degrees for 
limit cycle oscillations to appear. It was found that the frequency of 
the oscillations is Q = w; therefore, the Bode Plot for the simplified 
system would be as shown in Figure 23. 
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FIGURE 22. BLOCK DIAGRAM FOR PERTURBED r c C-SYSTEM 
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It would appear from Figure 23 that lead, lag, or lead-lag compensa¬ 
tions are not the answer to the problem since the phase shift will 
always go from -90 degrees to -270 degrees passing through -180 degrees. 
This was verified by analog simulation using different compensations. 

The amplitude and frequency of the limit cycles could be greatly 
altered but the oscillations, however small, still remained. Two fixes 
remained promising: 1) increase the deadband on the relay so that the 
oscillations present in the error signal would be within the relay 
deadband and would not drive the relay; and 2) alter the system so the 
phase shift never became greater than -180 degrees and thus limit cycles 
could not occur. 

Figure 13 shows the inherent limit cycle oscillations characteristic 
of the r q 0-simplified-controller. When the relay deadband was increased 
from e Q = .02 to e Q = .08, the limit cycles completely disappeared. The 
system became stable and transients resulting from a step elevator 
deflection soon died out. Figure 24 shows the r^e-system with e Q = .08 
and all other parameters remaining the same as in Figure 13; limit 
cycle oscillations are no longer present in the system. Although this 
fix works perfectly well, it is undesirable to make the oscillation free 
operation of the adaptive control system depend on the deadband size of 
the relay. If the deadband were made sufficiently large to eliminate 
limit cycles at the largest input command possible for safe aircraft 
operation, then the controller would be relatively nonresponsive at 
small input commands. This fine balance of deadband size would be most 
difficult to achieve and is therefore undesirable. 

The second fix, altering the basic system to avoid making the phase 
shift pass through -180 degrees seemed a reasonable approach. If the 
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second order servo were taken out of the perturbed r^C-control loop 
shown in Figure 22, then the phase shift would go from 0 degrees to 
-90 degrees instead of from -90 degrees to -270 degrees and the phase 
shift would never pass through -180 degrees. Thus, the system would 
not exhibit limit cycle oscillations. This approach is considered in 
the next chapter. 
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V. MODIFICATION OF THE SIMPLIFIED r,C-CONTROLLER 

0 


If a second nonvarying-C* control system could be developed similar 
to the first whose phase shift did not pass through -180 degrees, but 
whose variable gains approached their ideal values as the system error 
approached zero, then the new controller would be adaptive and would 
not exhibit limit cycles. Again, a simplified version of the complete 
system was studied first due to the complex nature of the system 
equations. The block diagram for the proposed altered r fi C-system is 
shown in Figure 25. The set of equations which describe this system 
are 


e = 

o e 

£ = Sr i - Ve = [6 s r { - <S { ]« e 


dr, 


at- = 


5 = r.c 

e o 


6 + 2^co6 + co 2 6 = co 2 <5 . 

e e e e 


Using the Gradient Technique, set 

2v “ It [ «6 r 6 ' 6 « ]! 


which comes from the error equation. Taking the derivative yields 
, , _ dr. 

dt = K7 ^ 6 6 r 6 " dt - ' 

dr 6 6 dr 6 

If is now chosen to be = K^6 e sgne, then the gradient becomes 


dt ^6 r 6 “ ^6^ e 6 <S e Sgn£ 

HF = h e sgne = ^ ,e| • 
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FIGURE 25. BLOCK DIAGRAM FOR MODIFIED f^C-SYSTEM 






































Since 8 g is always negative, the gradient is therefore also always 
negative and thus v and e approach zero as time proceeds. The ideal 
value of the variable gain may also be calculated from the error 
relationship: 


e = er. - 6.6 

o o e 

. . r.* = — for e s 0 . 

6 36 

The new simplified r^C-system is therefore adaptive since the variable 
gain will approach its ideal value as the error approaches zero and it 
seems to be a reasonable system to eliminate the parasitic oscillations 
of the previous T^C-controller. 

A short study of the stability of the modified r^C-controller will 
show whether or not limit cycle oscillations should be expected to 
occur in the system. The set of equations which describe the controller 
may be combined to yield 


dr 6 


K 5 6 e 2(6 6 r 6 


3 5 ) 


and 


6 + 2coj6 + oj 2 6 = w 2 r x C . 

e e e o 


If C is set equal to C Q and the steady state operation is observed, 

then f* = 0, 6 =0, and 6 =0. Therefore, 
o e e 

0 = Ve 2 ( VS 0 + 

and 6 = C r. . The steady state values of r~ and 6 may be calculated 

e 0 o e 

0 0 


to give 


6, 3* 

— and 5 = C . 

5 o 86 % 0 


Now, substituting the small perturbation relationships 
6=6 +6 


e e Q e 

/\ 

r = r + r 
*6 1 6 q 1 6 
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into the combined controller equations yields 

+ *.)*[B 4 (r 4 ♦? 4 >-b 4 ] 

o o 

A A 

6 + 2cw6 + oj 2 (6 + 6 ) = oj 2 C (l\ + fj. 

6 6 6 6 O 0 0 

o o 

If the previous steady state relationships are used, the equations 
reduce to 

dr S 

dt = Ve 6 6 r 6 

O 


6 + 2^0)6 + u 2 6 = w 2 r f 

e e e o 


which results in the following characteristic equation, 
(X - Kg<5 e 2 3g)(X 2 + 2twX + to 2 ) = 0. 


The system has three roots and since 6^ is always negative, it is 
obvious that all three roots have negative real parts, which indicates 
that the system is always stable. Any disturbance to the system will 
die out and limit cycles cannot occur. 

Observing the perturbed r.C-system shows that the phase shift can 
never cross -180 degrees which also indicates limit cycles cannot occur. 
These observations will be verified by analog simulation of the control¬ 
ler. The perturbed system equations are 


C = C o 



e = ® r 6 " * 

where y and e drive the system until y = 0, e = 0, and r g = r g *. It is 
also known that 

6 = r.C = J- C + yC 

e 6 0 3 o o 
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and 6 = 6 + <5 where 6 = ■=- C . A relationship for pitch accel- 

e e 0 e e o 0 

eration may be derived from the plant equation 

^<5 - /s - ^ 

9 = 6 S 6 e = B 6 C 0 TT + M- = + M 


3 


6 e o M 6 p 6 e 


and introduced into the error equation to yield 

(3 3 

6 - <r* Y)(c ° e « + 'A 1 ' 6 « ( r c ° + 

e = y3 6 C o + Y 
e = Y3 5 C 0 . 

Substituting the error relationship into the equation for the variable 
gain results in the following derivation: 


dr 6 d A 


dY 




& ■ K d y 1 r Vo 
°6 


and introducing the relay, 

& = (Y6 { C 0 ). 

3 <S 

This equation plus a relationship for the second order servo 

A ? 

6 + 2z;m<5 + oi 2 6 = w 2 yC describes the perturbed modified IAC-system 
shown in Figure 26. This figure shows that in the perturbed control 
loop, the maximum possible phase shift is -90 degrees. Therefore, the 
phase shift cannot pass through -180 degrees and limit cycles cannot 
occur. 

The r^C-controller was shown to exhibit limit cycles in Figures 7 
and 8. With the same parameters for the second order servo and variable 
gain, the modified r^C-controller is shown in Figure 27. The system 
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FIGURE 26. BLOCK DIAGRAM FOR PERTURBED MODIFIED r c C-SYSTEM 
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FIGURE 27. ANALOG TRACE FOR MODIFIED FgC-SYSTEM, NO LIMIT CYCLES 




























































































































































































































































does not exhibit limit cycle oscillations as the variable gain 
approaches its ideal value and the error approaches zero. A more 
detailed study of the control systems will be contained in Chapter VII 
for the complete nonvarying-C* adaptive control system. 
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VI. THE MODIFIED NONVARYING-C* ADAPTIVE CONTROLLER 


The simplified controller studied in the previous chapter may be 
expanded to the complete modified nonvarying-C* adaptive control 
system. The system must have the same fundamental aircraft equation 
as the adaptive control system studied in Chapter II, 

0=B9 + Bn+3p6 . 

q n o e 

The other governing equations for the controller were developed by 
trial and error using the information available from the simplified 
system and the Gradient Technique to insure that the system was 
adaptive. The error equation was computed for each configuration and 
from this relationship, the gradient could be set up. By choosing 
desirable values for the variable gains, the sign of the gradient 
could be calculated, and if the gradient was negative for all time, 
then the system error would approach zero. Using this method, the 
equations for the modified nonvarying-C* adaptive controller were 
calculated: 

E = r 6 [e + e 5 (r n n + r q e) - (e q e + e n n)] - 8 s « e 
dr, 

ar = V sgn E 

dr 

JT = K „ r « " sgn E 

dr 

= K r o0 sgn e 
dt q { 3 

<5 = r.[c + r n + r 0 ] 

e o n q 

6 + 2?(jo6 + o) 2 6 = (jo 2 6 

e e e e 
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The block diagram for this system is shown on Figure 28. Rearranging 
the error equation, 

6 ■ + r j te q ♦ s 6 r q - 6q ]8 

+ r A * 6 s r n - 8 n ]n , i 

allows the calculation of the ideal values of the variable gains: 


for e = 0 




b - e 

it H g 


0; 


B - B 

* _ n n 


n 


From the previous error relationship, also consider. 


dv 

dt 




1 


dr, 


dr. 


t r 6 6 6 - e s ] Hr + ^ + e 6 r q - 6 q ] df 

dr ) 


l 


*t K* Vn - ^ 

n 


n J dt 


Choosing 


dr. dr . dr 

i. - = K.<5 G ; Tj~ = K r.0G ; -rr^- = K r.nG 

dt 5 e dt q 6 dt n 5 

yields ^ = B^gG which is negative; therefore e approaches zero. From 
the fundamental aircraft equation and the relationship for 6 (neglecting 

S 

the servo) it may be shown that 

« ■ V + V + [r S (C + V + r n" )] ' 

If the ideal values of the variable gains are used in the previous 
equation. 


.. - . - b, B n - 3 n . 

6 = B n e + Bn + B, -J- [c + —3- 0 + 

q n 6 B 5 


B - B 
-As—- n] 

p* 
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FIGURE 28. gLOCK DIAGRAM FOR MODIFIED 
NONVARYING-C CONTROL SYSTEM 












































































then 0 = $q0 + 3 n n + 3^C, which is the C*-critenon as developed in 
Chapter II. The system has thus been proven to be adaptive and when 
the error is zero the aircraft behaves as if the fundamental aircraft 
equation were exactly that of the criterion. A detailed comparison of 
the two nonvarying-C* adaptive controllers will be shown in the next 
chapter. 
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VII. ANALOG STUDY OF TWO NONVARYING-C* ADAPTIVE CONTROL SYSTEMS 


The analog diagram used for the simulation of the original C*- 
controller was shown in Figure 14. The analog diagram for the modi¬ 
fied C*-controller is shown in Figure 29. The potentiometer settings 
for the aircraft and controller parameters for these analog diagrams 
are listed on Table I. All analog simulation was done on the Cl 5000 
analog computer. 

The original controller was studied in Chapter III for a flight 
condition of Mach 2 at 50,000 feet. Analog traces confirm the values 
of the variable gains listed on Table II without the complications 
introduced by the servo or the actuator. Traces are also shown for the 
control system with the actuator and with both the actuator and the 
servo in the system, resulting in limit cycles. Similar studies were 
made at flight conditions of Mach .9 at sea level and Mach .2 at sea 
level with comparable results. The flight parameters for the simula¬ 
tion were taken from F-4 data which appears on Table I. 

The modified controller was also studied at the same three flight 
conditions. The aircraft parameters were for the same aircraft and 
also appear on Table I. The ideal values for the variable gains of the 
modified controller are shown on Table II. The traces for the free 
aircraft without and with the actuator will be the same as shown in 
Figure 15 and 16 respectively since the fundamental aircraft equations 
are the same for both controllers. The modified controller, with 
actuator and feedback, for flight condition Mach 2 at 50,000 feet is 
shown in Figure 30. Note that the values of the variable gains approach 
their ideal values as listed in Table II even though the complication 
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FIGURE 29. ANALOG DIAGRAM FOR MODIFIED NONVARYING-C CONTROL SYSTEM 










































































FIGURE 29 . ANALOG DIAGRAM FOR MODIFIED NONVARYING-C CONTROL SYSTEM 
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FIGURE 29. ANALOG DIAGRAM FOR MODIFIED NONVARYING-C CONTROL SYSTEM 

















































































































































































































































































































































































introduced by the actuator is present. This results from the fact that 
the actuator, like the second order servo, is no longer in the perturbed 
system control loop and does not effect the values attained by the 
variable gains. Figures 31, 32, and 33 show the modified controller 
with the. actuator and the second order servo at flight conditions of 
Mach 2 at 50,000 feet, Mach .2 at sea level, and Mach .9 at sea level 
respectively. Limit cycle oscillations did not appear at any flight 
condition, further the controller was completely stable at all flight 
conditions studied. The modified nonvarying-C* adaptive control system 
was thus shown to be practical and free of limit cycle oscillations by 
analog simulation. 
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FIGURE 32. ANALOG TRACE FOR C -MODIFIED CONTROLLER, FLIGHT CONDITION 

MACH .2, SEA LEVEL C^.05 C 2 =.012 
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FIGURE 33. ANALOG TRACE 


FOR C -MODIFIED CONTROLLER, FLIGHT CONDITION 
MACH .9, SEA LEVEL C =.05 C 2 =.0l2 








































































































































































































































































VIII. CONCLUSIONS AND RECOMMENDATIONS 


The original nonvarying-C* adaptive controller had inherent limit 
cycle oscillations in the system. The two fixes for the problem were 
to increase the deadband on the relay and to modify the system to 
prevent the phase shift from passing through -180 degrees, both of 
which eliminated the oscillations. Modification of the controller, 
being the most practical approach, was studied thoroughly. 

A modified nonvarying-C* adaptive controller was designed and 
shown to be practical by analog simulation. An analytic study of the 
stability of the modified system together with a study of the perturbed 
control loop resulted in proof that limit cycles could not occur in the 
modified controller even under large input commands. The complications 
introduced by the actuator and the servo in the modified controller did 
not prevent the variable gains from approaching their ideal values as 
in the original controller, since the actuator and servo are not in the 
perturbed control loop of the modified system. 

Future study of the modified controller might include: an investi¬ 
gation of optimal gain constants, simplification of the modified system 
by further analog study, effects of gust inputs, a study of signal 
noise, a study of elevator hysteresis, and a control system flight 
simulation. The gains or constants K^, K^, and K n determine how fast 
or slow the integrations resulting in r g , r q , and r p take place. An 
optimal value for each of these constants will result in a rapid 
approach of the variable gains to their ideal values without a large 
amount of initial overshoot at all flight conditions. The quicker the 
variable gains approach their ideal values for a specific flight 
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condition, the quicker the system error becomes zero and the aircraft 
will be flying as if its governing equation were exactly that of the 
C*-criterion. 

The modified controller could probably be simplified from its present 

configuration. A more intensive analog study of the system would reveal 

whether or not the modified controller could be simplified and retain 

comparable adaptive characteristics. Figure 28 shows that is a 

common multiplication factor in the calculation of both r and r . 

q n 

By adjusting the constants and K n , this multiplication could proba¬ 
bly be eliminated. Other simplifications might result from continued 
analog simulation resulting in a more compact controller with acceptable 
characteristics. 

A complete study of signal noise in the system, elevator hysteresis, 
and the effects of gust inputs should be made on the modified control¬ 
ler. These complications are all present in a real physical control 
system and should be studied by simulation on the analog computer. 

A step gust input and signal noise in the mechanical linkages could be 
simulated and summed with the angle of attack signal to give total angle 
of attack; whereas, a model for elevator hysteresis could be patched in 
series with the actuator and series servo to give the elevator signal. 
Previous control systems have been shown to operate perfectly without 
the above conditions, but to fail as an operational adaptive controller 
with the introduction of one or more of these complications. 

After the system appears to be performing satisfactorily by the 
analog simulation, a flight simulation could be performed to determine 
pilot opinion of the controller. The Naval Postgraduate School variable 
stability flight simulator is capable of patching directly to the hybrid 
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computer's analog board. By so doing, the flying qualities of the air¬ 
craft with the modified controller could be studied at different flight 
conditions simply by resetting the potentiometers or by using a digital 
program to set the potentiometers for each flight condition or for each 
different aircraft under study. Since the Naval Postgraduate School 
has many naval aviators, the flight simulator could be flown by a 
number of pilots and rated to determine handling qualities of the air¬ 
craft with the modified controller. 
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